Flows and heterogeneities with a vane tool: MRI measurements 
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Synopsis 



We study the local flow properties of various materials in a vane-in-cup geometry. We 
use MRI techniques to measure velocities and particle concentrations in flowing Newtonian 
Q I fluid, yield stress fluid, and in a concentrated suspension of noncolloidal particles in a yield 

C/3 . stress fluid. In the Newtonian fluid, we observe that the 0-averaged strain rate component 

dre decreases as the inverse squared radius in the gap, in agreement with a Couette analogy. 
This allows direct comparison (without end-effect corrections) of the resistances to shear 
in vane and Couette geometries. Here, the mean shear stress in the vane-in-cup geometry 
'^ . is slightly lower than in a Couette cell of same dimensions, and a little higher than when 

^ ' the vane is embedded in an infinite medium. We also observe that the flow enters deeply 

O , the region between the blades, leading to signiflcant extensional flow. In the yield stress 

fluid, in contrast with the usually accepted picture based on simulation results from the 
literature, we find that the layer of material that is sheared near the blades at low velocity 
\^ ' is not cylindrical. There is thus a significant extensional component of shear that should be 

taken into account in the analysis. Finally and surprisingly, in the suspension, we observe 
,—1. I that a thin non-cylindrical slip layer made of the pure interstitial yield stress fluid appears 

^^ , quickly at the interface between the sheared material and the material that moves as a rigid 

lO ' body between the blades. This feature can be attributed to the non-symmetric trajectories 

f»^ . of the noncolloidal particles around the edges of the blades. This new important observation 

f^ ' is in sharp contradiction with the common belief that the vane tool prevents slippage, and 

^D , may preclude the use of the vane tool for studying the flows of pasty materials with large 

particles. 

rS : I Introduction 



Experimental investigations of the rheology of concentrated suspensions often involve a vane- 
in-cup geometry (see Barnes and Nguyen (2001) for a review). The vane tool offers two main 
advantages over other geometries. First, it allows the study of the properties of structured 
materials with minimal disturbance of the material structure during the insertion of the tool 
[Dzuy and Boger (1983); Alderman et al. (1991)]. It is thus widely used to study the properties 
of gels and thixotropic materials [Alderman et al. (1991); Stokes and Telford (2004)] and for in 
situ study of materials as e.g. in the context of soil mechanics [Richards (1988)]. Second, it is 
supposed to avoid wall slip [Keentok (1982); Dzuy and Boger (1983); Saak et al. (2001)], which 
is a critical feature in concentrated suspensions [Coussot (2005)]; the reason for this belief is 
that the material sheared in the gap of the geometry is sheared by the (same) material that 
is trapped between the blades. Consequently, it is widely used to study the behavior of pasty 
materials containing large particles, such as fresh concrete [Koehler et al. (2006); Estelle et 
al. (2008); Wallevik (2008); Jau and Yang (2010)] and foodstuff [Stokes and Telford (2004); 
Martmez-Padilla and Rivera- Vargas (2006)]. 
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The constitutive law of materials can be obtained from a rheological study with the vane- 
in-cup geometry provided one knows the coefficients - called "geometry factors" ~ that allow 
the conversion of the raw macroscopic data (torque, rotational angle or velocity) into local data 
(shear stress, shear strain or shear rate). However, in contrast with other classical geometries, 
even the a priori simple linear problem (for Hookean or Newtonian materials) is complex to solve 
with a vane tool. This linear problem was studied theoretically by Sherwood and Meeten (1991) 
and Atkinson and Sherwood (1992) in the general case of a A^-bladed vane tool embedded in an 
infinite linear medium. The analytical expression found for the torque vs. rotational velocity is 
in rather good agreement with macroscopic experimental data [Sherwood and Meeten (1991)]. 
Note however two possible shortcomings of this theoretical approach for its use in practice: the 
blades are infinitely thin and there is no external cylinder. 

There is no such approach in the case of nonlinear media (i.e. complex fluids). A practical 
method used to study the flow properties of non-linear materials, known as the Couette analogy 
[Bousmina et al. (1999); Ait-Kadi et al. (2002); Estelle et al. (2008)], consists in calibrating the 
geometry factors with Hookean or Newtonian materials. One defines the equivalent inner radius 
Ri^eq of the vane-in-cup geometry as the radius of the inner cylinder of a Couette geometry that 
would have the same geometry factors for a linear material. For any material, all macroscopic 
data are then analyzed as if the material was sheared in a Couette geometry of inner cylinder 
radius Ri^cq. The nonlinearity (that affects the flow field) is sometimes accounted for as it is in a 
standard Couette geometry [Estelle et al. (2008)]. This approach may finally provide constitutive 
law measurements within a good approximation [Baravian et al. (2002)]. 

However, simulations and observations show that Ri^eq is not a universal parameter of the 
vane tool independent of the properties of the studied material. While the streamlines go into the 
virtual cylinder delimited by the blades in the case of Newtonian media [Baravian et al. (2002)], 
yielding an equivalent radius lower than the vane radius [Sherwood and Meeten (1991); Atkinson 
and Sherwood (1992)], it was found from simulations [Barnes and Carnali (1990); Savarmand 
et al. (2007)] that the streamlines are nearly cylindrical everywhere for shear-thinning fluids if 
their index n is of order 0.5 or less, and thus that Ri^eq = Ri in these cases. Moreover, for 
yield stress fluids, simulations and photographs of the shearing zone around a four-bladed vane 
rotating in Bingham fluids [Keentok et al. (1985)], simulations of Herschel-Bulkley and Casson 
fluids flows in a four-bladed vane-in-cup geometry [Yan and James (1997)], and simulations of 
Bingham fluids flows in a six-bladed vane-in-cup geometry [Savarmand et al. (2007)], all show 
that at yield (i.e. at low shear rates), the material contained in the virtual cylinder delimited 
by the blades rotates as a rigid body, and that it flows uniformly in a thin cylindrical layer 
near the blades. This is now widely accepted [Barnes and Nguyen (2001)] and used to perform 
a Couette analogy with Ri^cq = Ri', the yield stress Ty is then simply extracted from torque 
T measurements at low velocity thanks to Ty = T / {27: HRf), where H is the vane tool height 
(neglecting end effects) [Nguyen and Boger (1992)]. 

The flow field in a vane-in-cup geometry and its consequences on the geometry factors have 
thus led to many studies. However, only theoretical calculations, macroscopic measurements 
and simulation data exist in the literature: there are no experimental local measurements of the 
flow properties of Newtonian and non-Newtonian materials induced by a vane tool except the 
qualitative visualization of streamlines made by Baravian et al. (2002) for Newtonian media, 
and the photographs of Keentok et al. (1985) for yield stress fluids. Moreover, while the main 
advantage of the vane tool is the postulated absence of wall slip, as far as we know, this widely 
accepted hypothesis has been neither investigated in depth nor criticized. In order to provide 
such local data, we have performed velocity measurements during the flows of a Newtonian 
medium and of a yield stress fluid in both a coaxial cylinder geometry and a vane-in-cup geome- 
try. We have also performed particle concentration measurements in a concentrated suspension 
of noncolloidal particles in a yield stress fluid, which is a good model system for complex pastes 
such as fresh concrete [Mahaut et al. (2008a, b)]. Our main results are that: 



(i) in the Newtonian fluid, the 0-averaged strain rate component dre decreases as the inverse 
squared radius in the gap, as in a Couette geometry, which aUows direct determination 
(without end-effect corrections) of the value of Ri^eq'- it is here found to be lower than Ri, 
but slightly higher than for a vane in an infinite medium; the flow enters deeply the region 
between the blades, leading to a significant extensional flow; 

(ii) in the yield stress fluid, in contrast with results from the literature, the layer of material 
that is sheared near the blades at low velocity does not have a cylindrical shape; 

(iii) in the suspension of noncolloidal particles in a yield stress fluid, the noncolloidal particles 
are quickly expelled from a thin zone near the blades, leading to the development of a thin 
slip layer made of the pure interstitial yield stress fluid, in sharp contradiction with the 
common belief that the vane tool prevents slippage. 

In Sec. II, we present the materials employed and the experimental setup. We present the 
experimental results in Sec. Ill: velocity profiles obtained with a Newtonian oil and with a yield 
stress fluid are presented in Sec. IIIA and Sec. IIIB, while Sec. IIIC is devoted to the case of 
suspensions, with a focus on the slip layer created by a shear-induced migration phenomenon 
specific to the vane tool. 

Throughout this paper, we use cylindrical coordinates {r,9,z). All fiows are supposed to be 
z invariants (i.e. there are no flow instabilities). We deflne the ^-average /(r) of a function 
/(r,0)as/(r) = (l/27r)/(fV(r,e)d0. 

II Materials and methods 

A Materials 

We study three materials: a Newtonian fluid, a yield stress fluid, and a concentrated suspension 

of noncolloidal particles in this yield stress fluid. 

The Newtonian fluid is a silicon oil of 20 mPa.s viscosity. 

The yield stress fluid is a concentrated water in oil emulsion. The continuous phase is dodecane 

oil in which Span 80 emulsifier is dispersed at a 7% concentration. A 100 g/1 CaCl2 solution is 

then dispersed in the oil phase at 6000 rpm during 1 hour with a Sliverson L4RT mixer. The 

droplets have a size of order 1 ^m from microscope observations. The droplet concentration is 

75%, and the emulsion density is pf = 1.01 gcm~^. The emulsion behavior, measured through 

coupled rheological and MRI techniques described in Ovarlez et al. (2008) (see Fig. 1), is well 

fltted to a Herschel-Bulkley behavior t = Ty + rjHBi"' with yield stress Ty = 20.6 Pa, consistency 

'nHB=Q-8 Pas'^'^^, and index n = 0.44. 

The suspension is a suspension of monodisperse polystyrene beads (density pp = 1.05 gcm~^, 
diameter d = 250 ^um) suspended at a 40% volume fraction in the concentrated emulsion de- 
scribed above. The density matching between the particles and the yield stress fluid is sufficient 
to prevent shear-induced sedimentation of the particles in the yield stress fiuid [Ovarlez et al. 
(2010)]; in all experiments, we check that the material remains homogeneous in the vertical 
direction by means of MRI density measurements. 

B Rheometry 

The rheometric experiments are mainly performed within a six-bladed vane-in-cup geometry 
(vane radius Ri = 4.02 cm, outer cylinder radius Rq = 6 cm, height H =11 cm). The shaft 
radius is 1.1 cm and the blade thickness is 6 mm. Other experiments are performed with a 
wide-gap Couette geometry of slightly different inner cylinder radius Ri = 4.15 cm, due to the 
presence of sandpaper (the other dimensions were identical). The inner cylinder of the Couette 
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Figure 1: Constitutive law of the emulsion measured locally in the gap of a Couette cell through 
MRI techniques [Ovarlez et al. (2008)]. The empty squares are local data; the solid line is a 



Herschel-Bulkley fit to the data r 
n = 0.44. 



Ty + rjHBJ"' with Ty = 20.6 Pa, riHB=G-8 Pas' 
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geometry, and the outer cylinder of both geometries are covered with sandpaper of roughness 
equivalent to the size of the largest elements of the materials studied in order to avoid wall slip. 
In the rheometric experiments presented here, we control the rotational velocity of the inner 
cylinder, with values ranging from 0.1 to 100 rpni. 



C MRI 

Proton MRI [Callaghan (1991)] was chosen as a non-intrusive technique in order to get measure- 
ments of the local velocity and of the local bead concentration inside the sample. Experiments 
are performed on a Bruker 24/80 DBX spectrometer equipped with a 0.5T vertical supercon- 
ductive magnet with 40 cm bore diameter and operating at 21 MHz (proton frequency). We 
perform our experiments with a home made NMR-compliant rheometer, equipped with the ge- 
ometries described in the previous section. This device was already used in a number of previous 
rheo-nmr studies [Raynaud et al. (2002); Rodts et al. (2004); Ovarlez et al. (2006)], and is fully 
described in Raynaud et al. (2002). The volume imaged is a (virtual) rectangular portion of 
4 cm in the axial (vertical) direction with a width (in the tangential direction) of 1 cm and a 
length of 7 cm (in the radial direction, starting from the central axis). Velocity and concentra- 
tion profiles Vg{r) and 4>{r), averaged over the vertical and tangential directions in this volume, 
are obtained with a resolution of 270 nm in the radial direction. This volume is situated at the 
magnet center (so as to minimize the effects of field heterogeneities) and sufficiently far from 
the bottom and the free surface of the rheometer so that ffow perturbations due to edge effects 
are negligible. We checked that the velocity and concentration profiles are homogeneous along 
the vertical direction in this volume, which justifies averaging data over this direction. 

Details on the sequence used to obtain velocity profiles can be found in [Raynaud et al. 
(2002); Rodts et al. (2004)]. While it is possible to get 2D or 3D maps of 2D or 3D velocity 
vectors [Rodts et al. (2004)], such measurements may actually take minutes and require complex 
synchronization of the MRI sequences and of the geometry position. However, we will show in 
the following that the azimuthal velocity alone provides a valuable information that can be 
sufficient for most analyses; in particular it allows computation of the ^-averaged strain rate 
component dre and thus the ^-averaged shear stress fro (and the torque T) in the case of the 
Newtonian oil. That is why we have chosen to limit ourselves to ID profiles of ID velocity 
measurements, namely the azimuthal velocity V0{r,t) as a function of the radius r and time 



t, for which a single measurement may take as httle as 1 s; this has allowed us to perform a 
sufficient number of experiments, with various materials, geometries, and rotational velocities. 
Depending on the time over which this measurement is averaged as compared to 2tt/{NQ) 
where iV = 6 is the number of blades and Q is the rotational velocity, this measurement may 
provide either a time (or 0)-averaged azimuthal velocity Vg{r) = (l/lir) L^ Vg{r,9)d9 or an 
instantaneous (transient) azimuthal velocity Vo{r,t). In this latter case, the 9 dependence of 
the azimuthal velocity Vg{r,9) at a given radius r can then be easily reconstructed by simply 
replacing the time t dependence by an angular 9 dependence with 9 = Vtt. It should also be noted 
that due to incompressibility of the materials we study, the Vr{r,9) field can be reconstructed 
thanks to {l/r)dr{rVr) + {l/r)de{Ve) = with Vr{Ro,9) = 0; however, this derivation from the 
experimentally measured values of Vg{r,9) cannot be very accurate. Finally, from Vr{r,9) and 
Vg{r,9), we are also able to evaluate the strain rate components drr{r,9) = —d0g{r,9) = drVr 
and drg{r,9) = (1/2) [(l/r)3e(T4) +r3r(Ve/r)]. Note that the derivative dxf with respect to 
coordinate x of experimental data /(xj) measured at regularly spaced positions Xi was computed 
as: dxf{xi) = [/(xj+i) - f{xi-i)]/[xi+i - Xj_i]. 

The NMR sequence used in this work to measure the local bead concentration is a modified 
version of the sequence aiming at measuring velocity profiles along one diameter in Couette 
geometry [Hanlon et al. (1998); Raynaud et al. (2002)], and is described in full detail in Ovarlez 
et al. (2006). The basic idea is that during measurements, only NMR signal originating from 
those hydrogen nuclei belonging to the liquid phase of the sample (i.e. both the oil and water 
phase of the emulsions) is recorded: the local NMR signal that is measured is thus proportional 
to 1 — 0, where (j) is the local particle volume fraction. A rather low absolute uncertainty of 
±0.2% on the concentration measurements values was estimated in Ovarlez et al. (2006). All 
volume fraction profiles are measured at rest, after a given flow history. This is possible because 
the particles do not settle in the yield stress fluid at rest: the volume fraction profile induced 
by shear is gelled by the interstitial yield stress fluid. 

Ill Experimental results 

In this section, we study successively the flow properties of the Newtonian oil, the yield stress 
fluid, and the suspension. 

A Newtonian fluids 

In this section, we study the flows observed with a Newtonian fluid. We first present a basic 
theoretical analysis of the fiows in a vane-in-cup geometry as compared to flows in a standard 
Couette geometry, which provides the basis for a Couette analogy. The ^-averaged azimuthal 
proflles Ve{r) are then shown, and are compared to predictions of the Couette analogy. The full 
velocity fleld V0{r,9), Vr{r,9) is finally presented and analyzed. 

1 Couette analogy: theoretical analysis 

The stress balance equation projected along the azimuthal axis is: 

il/r)drir^Tre) + deiree) - dep = (1) 

where Tij is the deviatoric stress tensor and p the pressure. 
The strain rate tensor component d^e is given by: 

dre(.r,9) = ^[{l/r)de{Vr) + rdr{Ve/r)^ (2) 

We recall that the constitutive law of a Newtonian fluid of viscosity rj is: 

Tij = 2r]dij (3) 



In all the following analysis, we assume a no-slip boundary condition at the walls of the inner 
tool and of the cup. 

Couette geometry 

In a standard - coaxial cylinders - Couette geometry, due to cylindrical symmetry, Eq. 1 becomes 
dr{r'^Tj-.g) = which means that the whole shear stress distribution T^o{r) in the gap is known 
whatever the constitutive law of the material is. If a torque T{Vl) is exerted on the inner cylinder 
driven at a rotational velocity O, Tro^r) is given by: 

rMr) = -^ (4) 

For a Newtonian fluid of viscosity ij, it follows that the strain rate component droir) is given by: 



As Eq. 2 becomes dreir) = (1/2) r driVe/r) with cylindrical symmetry, due to the boundary 
conditions Ve(i?j) = QRi and Vq^Rq) = 0, from fj^^2dr0{r)/rdr = i7, one gets alternatively 
dre{^) = ~^^f-Ro/[^^(-^o ~ ^f)]- This yields the following azimuthal velocity profile: 

p2 p2 _ 2 

Finally, the viscosity 77 of a Newtonian fluid is obtained from the measured torque/rotational 
velocity relationship T{Q) through 

^ n attHrir] ^^ 

These equations will be used for the comparison with the flows observed in a vane-in-cup 
geometry, in particular to determine the radius Ri^cq of the equivalent Couette geometry. 

Vane-in-cup geometry 

In a vane-in-cup geometry, there is a priori no cylindrical symmetry and all quantities a priori 
depend on 9. However, averaging Eq. 1 over 6 yields dr{r'^fre) = 0. This means that what 
is true in a Couette geometry, Tre{r) = Tre{Ri)R^ /r"^, is still true on average with a vane-in- 
cup geometry: fre{f) = fre{Ri)Rf/r'^ independently of the material's constitutive law. Note 
that this derivation is true only between Ri and Rq', this is not true for the material between 
the blades as the unknown Tij distribution in the blades contributes to the ^-average. The 
link between this stress distribution and the torque T{Q) exerted on the vane tool may then 
seem difficult to build. However, it can be equivalently computed on the outer cylinder as 
T = j'^"" Tre{Ro,0)HRlde = 2tt HRlfre{Ro)- This means that Eq. 4 is still valid for the 9- 
averaged shear stress in the vane-in-cup geometry, for Ri < r < Rq'. 

^Mr) = -^ (8) 

Prom the ^-averaged Eq. 3, this means that the ^-averaged strain rate distribution in a 
Newtonian fluid, for Ri < r < Rq, is: 



From the 0-averaged Eq. 2 dro = {l/2)r dr{Ve/r), this means that the ^-averaged azimuthal 
velocity profile of a Newtonian fluid of viscosity r/ in a vane-in-cup geometry for Ri < r < Rq, 
with a boundary condition VelRo) = is given by: 

Voir) = ^M^i^ (10) 

Finally, the only difference with a standard Couette flow, as regards these ^-averaged quan- 
tities, is that we do not know the value of Vo{Ri); we only know that Ve(i?i, 2TTk/n) = QRi, for 
k integer, where A^ is the number of blades. This means that dre{r) and VqIv) follow the same 
scaling with r and 0, as in the standard Couette geometry, but with a different prefactor. 

Nevertheless, these equations provide a new insight in the Couette analogy. The usual way 
of performing the Couette analogy consists in defining the radius of the equivalent Couette 
geometry Ri^eq as the radius that allows measuring the viscosity r/ of a Newtonian fluid with 
the standard Couette formula. From Eq. 7, ij should then be correctly obtained from the 
torque/rotational velocity relationship T{Q) measured in a vane-in-cup geometry with: 

_ T{n) Rl - Rleg 
"- n A.HRlRl^ ^''^ 

From Eqs. 7 and 11, it means in particular that the torque Tyanc exerted by the vane tool is 
decreased by a factor 

T R^ 1 _ p2 / p2 



^Couette ^j 1 - ^i,eq/^o 

as compared to the torque Tcouette exerted by the inner cylinder of a Couette geometry of same 
radius Ri at a same rotational velocity. 

Here, from Eqs. 10 and 6, we see that from the local flow perspective, there is a Couette 
analogy in the sense that the 0-averaged azimuthal velocity (and shear) profiles will be exactly 
the same as in a Couette geometry. This defines a radius Ri^cq of the equivalent Couette geometry, 
such that Ve(r) and dre{r) for Ri < r < Rq are given by: 

Mr) = ^^#^ (13) 

jIq ^,eq 

Of course, these two definitions of Ri^cq are equivalent: combining Eqs. 10 and 11 yields Eq. 13. 
This point of view provides an additional meaning to the Couette analogy, namely the 
similarity of the average flows, and offers a new experimental mean to determine Ri^cq, which is 
more accurate than calibration. In rheological measurements, the T{i^) relationship has to be 
corrected for end effects [Sherwood and Meeten (1991); Martinez-Padilla and Quemada (2007); 
Savarmand et al. (2007)] and the Couette analogy has to be calibrated on a reference material 
of known viscosity. Here, the V0{r) or droir) measurements provide the value of Ri^cq directly 
without any correction, as only shear in the gap is involved in the analysis, and independent of 
the viscosity of the material. This will be illustrated in the following. 

Vane in a finite cup vs. vane in an infinite medium 

The only theoretical prediction of the stress field associated with a vane tool is that of Atkinson 
and Sherwood (1992) for an infinite A'^-bladed vane embedded in an infinite linear medium. 



In this case, it is shown that the torque T^anci^,N,Ri,Ro = oo) exerted on the vane is well 
approximated by 

T^!,nci^,N,Ri,Ro = Oo) 1 



?Coucttc {^, Ri,Ro = oo) N 



where Tcoucttc i^i Ri, Ro = oo) is the torque exerted on a cylinder of same radius Ri as the vane 
in an infinite medium (i.e. with Rq = oo). Eq. 15 is in agreement with experimental results 
[Sherwood and Meeten (1991)]. 

Sherwood and Meeten (1991) argue that, as the stress distribution varies as 1/r^ in a Couette 
geometry, {Ri/Ro)"^ should be of the order of 1% or less in order to nullify the influence of the 
outer boundary; this is clearly the case in their experiments and in field experiments where the 
vane is embedded e.g. in a soil; this is clearly not the case in our experiments and in most 
rheological experiments that make use of a vane-in-cup geometry. However, when the cup to 
vane radius ratio Rq/Ri is not large, no generic theoretical expression exists in the literature. 

Nevertheless, bounds of the value of the torque T^ane{^,N,Ri,Ro) can be derived using 
classical results of linear elasticity [Salengon (2001)]. Our starting points are the variational 
approaches to the solution of the Stokes equations describing the flow of an incompressible 
linear material induced by the rotation of an inner tool (of any shape) at a rotational velocity 
0, within a cup. In this framework, it can be shown that [Salencon (2001)]: 

/ 2ro-g„ dS - ^^ / T-jT-j dw < Tvane (^, N, Ri, Ro) < "^ / ^j^ij dw (16) 

where S^ denotes the inner tool- fluid interface and ujf the domain occupied by the fluid. In 
the first inequality, o"^ • is any stress tensor complying with the stress balance equations, a'^^ is 
the azimuthal component of the surface forces applied by the tool on the fluid and t-, is the 
deviatoric stress tensor associated to a'-. In the second inequality, d^- is the strain rate tensor 
associated with any velocity field V_' complying with the incompressibility constraint and the 
boundary conditions prescribed on the tool-fluid and cup-fluid interfaces. 
Eq. 16 leads in particular to the expected inequalities: 

Tvanc (^, N, Ri,Ro = Oo) < Tvanc i^, N, R4, Ro) < Tcoucttc i^, Ri, Ro) (17) 

The lower bound is obtained by using the velocity field defined by V_' = V_q n r r ^o^ r < Rq 
and V_' = for r > Ro, where V^qatr. r is the solution for the A^-bladed vane of radius Ri 
in a cup of radius Ro- V^ complies with the boundary conditions for the A^-bladed vane of 
radius Ri in an infinite domain problem. Then, putting this test velocity field within the second 
inequality (16) with Rq = oo and using 

T^anc i^, N, Ri,Ro) = 5 / d^jd,j doj (18) 

'^ Jw(N,Ri,Ro) 

yields the lower bound of the inequality (17) for the quantity T^iaic{^^N,Ri,Ro). In Eq. 18, 
dij denotes the strain rate tensor associated with V_q n r r while oj [N, Ri, Ro) is the domain 
occupied by the fluid. 

The upper bound is obtained using the test velocity field defined by V^ = Vg (r) eg with Vq 
defined by Eq. 6 for Ri < r < Ro and by Vg = 0,r for r < Ri. It is easily checked that V^ 
complies with the velocity boundary conditions for any vane-in-cup geometry with vane radius 
Ri and cup radius Ro- Putting this test velocity field into the second inequality (16) then yields 
the upper bound of inequality (17). 

Finally, combining inequalities (17), Eq. 15 and Eq. 7 yields 

1 _ iVl - :^^ < T..no {n,N,R,,Ro) ^ 

Njy RlJ - Tcouotte{n,R„Ro) - ^ ^ 



Of course, it is possible to improve the lower bound by determining admissible test stress 
fields for the problem under consideration and the inequalities (16). For example, let us consider 
the stress field defined between the two blades positioned at 6 = -im/N by 

f f \ f V \ Tfi ~\~ '2i I T \ 



with m> —1 for r < Ri and by 



2 

RiJ '""" ~"""' '"m + l\Ri 



a'^Q = -t[-^] ■ a00 = O ; a'„ = -t9——^ ( -^ ] (21) 



for Ri < r < Rq. This stress field complies with the balance equations within the fluid domain. 
Let us recall that a stress field does not need to be continuous to comply with the balance 
equations (of course, this stress field is not the solution of the problem). Putting this stress 
field into the first inequality (16) and using a numerical optimization tool to choose the optimal 
value of the parameter m yields a new lower bound for the A^-bladed vane in cup problem, 
which depends on N and Ri/Ro- In some cases, this test stress field improves the lower bound 
of inequality (19): e.g., for the geometry we use in this study (A^ = 6,Ro/Ri = 1-49), the new 
lower bound is 0.57 while the lower bound given by inequality (19) is 0.45. Nevertheless, such an 
improvement is not obtained for all parameter sets {N, Ri/Ro). It is thus necessary to compute 
the two lower bounds for each value of {N, Ri/Ro) in order to obtain the more accurate lower 
estimate of the torque. Although application of variational approaches to the derivation of 
estimates of the applied torque of a vane-in-cup problem is not classical, it is believed that 
such a strategy is able to provide useful results when no theoretical prediction of the solution is 
available for particular geometries. Lower bounds of Tvanc/^Coucttc computed using the approach 
presented above are displayed in Tab. 1 and are compared below to our results and to data in 
the literature. 

2 0-averaged profiles 

We first study the ^-averaged azimuthal velocity profiles Vg{r) observed during the flows of 
a Newtonian oil (Fig. 2). As shown above, these profiles can be used to check the validity 
of the Couette analogy and to determine the Couette equivalent radius Ri^eq- The azimuthal 
dependence of the velocity profiles between two adjacent blades of the vane tool will then be 
considered. 

In Fig. 2a we observe that the velocity profiles in the gap of a Couette geometry are, as ex- 
pected, in perfect agreement with the theory for a Newtonian fiow (Eq. 6). This first observation 
can be seen as a validation of the measurement technique. 

In the vane-in-cup geometry (Fig. 2b), we first note that the 0-averaged dimensionless az- 
imuthal velocity profiles Vg{r,0,)/0,Ri measured for several rotational velocities O are super- 
posed, as expected from the linear behavior of the material. We also remark that the material 
between the blades rotates as a rigid body only up to r ~ 3.1 cm, indicating that the shear 
fiow enters deeply the region between the blades (the vane radius is 4.02 cm). The whole 
limit between the sheared and the unsheared material in the (r, 6) plane will be determined in 
Sec. IIIA3 (Fig. 4b). We finally observe that the theoretical velocity profile for a Newtonian 
fiuid in a Couette geometry of radius equal to that of the vane lies above the data, as expected 
from the literature. This is also consistent with the observation that the shear fiow enters the 
region between the vane blades. 

In order to test the Couette analogy, we have chosen to plot the ^-averaged strain rate dro 
vs. the radius r in Fig. 3. This allows us to distinguish more clearly the difference between 
the experimental and theoretical flow properties than would the velocity profiles, because the 
velocity profile always tends to the same limit {V{Ro) = 0) at the outer cylinder whereas the 
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Figure 2: a) Dimensionless velocity profile Vg{r, Q)/nRi of a Newtonian oil in a Couette geometry 
{Ri = 4.15 cm), at various rotational velocities O ranging from 2 to 20 rpm; the solid line 
is the theoretical profile for a Newtonian fluid, b) Dimensionless ^-averaged velocity profile 
V0{r,Q)/VLRi of a Newtonian oil in a six-bladed vane-in-cup geometry [Ri = 4.02 cm) for i7 
ranging from 1 to 9 rpm; the vertical dashed line shows the radius of the vane; the dotted line 
is the theoretical profile for a rigid body rotation (for r < Ri); the solid lines are the theoretical 
profiles for a Newtonian fluid in Couette geometries of radii, from right to left: (i) Ri = 4.02 cm, 
(ii) Ri^cq = 3.90 cm, and (iii) Ri^th = 3.82 cm corresponding to the Atkinson and Sherwood 
(1992) theory in an infinite medium. 
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Figure 3: ^-averaged strain rate d^e vs. radius r for a Newtonian oil sheared at 1 rpm in a 
six-bladed vane-in-cup geometry. The vertical dashed line shows the radius of the vane. The 
solid lines are the theoretical strain rate profiles for a Newtonian fluid in Couette geometries of 
radii: (i) Ri = 4.02 cm (light grey), (ii) Ri^eq = 3.90 cm (black), and (iii) Ri^th = 3.82 cm (dark 
grey) corresponding to the Atkinson and Sherwood (1992) theory in a inflnite medium. 



strain rate profile does not. Note that velocity measurements could not be performed close to 
the blades, which explains why strain rate data are missing from 4 to 4.2 cm. 

In Fig. 3, we first note that dre is zero up to ~ 3 cm, which corresponds to the limit of 
the rigid motion of the material; dre then increases when r tends towards Ri as the material is 
more and more sheared between the blades. In the gap of the geometry, dro decreases when r 
increases. As expected, the theoretical strain rate profile for a Newtonian material in a Couette 
geometry of radius equal to that of the vane falls well above the data at any radius r (this 
was less obvious on the velocity profiles). We then observe that the data are well fitted to the 
theoretical strain rate profile (Eq. 14) for a Newtonian material flowing in an equivalent Couette 
geometry of inner cylinder radius Ri^eq = 3.90 cm {Ri^cq = 3.905 it 0.005 cm was obtained from a 
fit of the velocity profile to Eq. 13). This confirms that the ^-averaged strain rate dro decreases 
as the inverse squared radius in the gap, in agreement with the Couette analogy. 

From Eq. 12, we find Tvanc/Tcoucttc = 0.90 (let us recall that we do not need to consider end 
effects here because we determine the shear rate within the gap, and hence only the contribution 
to the torque from the material sheared in the gap). This value can now be compared to data 
from the literature. For a six-bladed vane tool in an infinite medium, the Atkinson and Sherwood 
(1992) theory would imply a theoretical Tvanc /^Couette = 0.83, which is 8% lower than what we 
measure, and corresponds to a theoretical "equivalent radius" Ri^th = 3.82 cm when the vane 
is embedded in a cup of radius Rq = 6 cm. Figs. 2b and 3 show that the flow characteristics 
predicted with this value of the equivalent radius can be distinguished from our experimental 
data and fall slightly below the data (a discrepancy could be expected as {Ri/Ro)'^ is not small 
in our experiment). 
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Table 1: Ratio Tvanc/Tcoucttc between the torque measured when straining a linear medium 
(viscous or elastic) in a vane-in-cup geometry and that measured in a coaxial cylinder geometry 
of similar dimensions, obtained in various theoretical, numerical and experimental studies of 
the literature; only data corrected for (or free from) end effects are shown. The number A^ of 
blades, the cup to vane radius ratio Ro/Ri, the shaft radius to vane radius ratio r^/Ri, and the 
blade thickness to vane radius ratio e/Ri, are displayed when provided in the manuscripts. The 
theoretical lower bound computed using variational approaches in Sec. IIIAl is also provided. 



We have gathered experimental and numerical data from the literature where the cup to vane 
radius ratio Ro/Ri is not large in Tab. 1; only data corrected for (or free from) end effects are 
shown. First, it should be noted that all torque data obey the theoretical inequalities computed 
using variational approaches in Sec. IIIAl. However, no clear trends emerge from the comparison 
of the data. The relative impact of the various geometrical parameters that may affect the flow 
field, namely the cup to vane radius ratio Ro/Ri, the shaft radius to vane radius ratio Vs/Ri, 
and the blade thickness to vane radius ratio e/i?j, cannot be determined at this stage. For 
example, in very similar geometries, Zhu et al. (2010) find a torque ratio Tyanc /^Coucttc close to 



that of Atkinson and Sherwood (1992) whereas Barnes and Carnah (1990) find a much lower 
torque ratio. The only noticeable difference between these two studies (apart from the numerical 
method) is that Vs/Ri is higher in Barnes and Carnali (1990), but our data, with a rather large 
value of Ts/Ri-, show different features. We actually note that our study is the only one to report 
a torque ratio higher than in an infinite medium; all other data report torque ratios up to 19% 
lower than expected in an infinite medium. In the general case of a finite vane-in-cup geometry, 
it thus seems that numerical investigations are still needed, and that, at this stage, a calibration 
has to be performed to get the geometry factors. We also expect that the bounds obtained using 
variational approaches in Sec. IIIAl can be improved. 

3 9 dependent profiles 

To better characterize the flow field, we now study the dependence of the velocity profiles on 
the angular position 6. We have performed experiments in which we measure one azimuthal 
velocity profile per second while the vane tool is rotated at 1 rpm, yielding 10 profiles between 
two adjacent blades. 



0.004 




0.000 



(a) 



120 



2 3 4 

Radius (cm) 



Figure 4: a) Azimuthal velocity profile V0{r,6) of a Newtonian oil sheared at 1 rpm in a six- 
bladed vane-in-cup geometry, for various angular positions, ^, between one blade (^ = ° ) and 
midway between adjacent blades {9 = 30 ° ). The vertical dashed line shows the radius of the 
vane. The dotted line is the profile for a rigid body rotation (for r < Ri) and the theoretical 
profile for a Newtonian fiuid in a Couette geometry of radius Ri (for r > Ri). b) Two-dimensional 
plot of the limit between rigid motion and shear (empty circles) for a Newtonian material in the 
six-bladed vane-in-cup geometry; the grey rectangles correspond to the blades. 



In Fig. 4a, we plot the velocity profiles Vg{r,9) measured at different angles 9. We first 
observe that the velocity profile which starts near from a blade tip (corresponding to ^ = ° by 
definition) is very different from the velocity profile in a Couette geometry of same radius: it 
starts with a much steeper slope, which means that the blades tip neighborhoods are regions 
of important shear as already observed by Barnes and Carnali (1990). We then observe that, 
as expected from the 0-averaged velocity profiles, the shear fiow enters more and more deeply 
the region between the blades as 9 tends towards 30 ° (corresponding to midway between two 
adjacent blades); at this angular position, the rigid rotation stops at Ri ~ 3.05 cm. From all the 
velocity profiles, we finally extract a 2D map of the limit Ri{9) between rigid rotation and shear, 
which is depicted in Fig. 4b. This provides an idea of the deviation from cylindrical symmetry, 
and will be compared in the following to the case of yield stress fluids. Note that eddies are likely 
to be present in the "rigid" region [Moffatt (1964); Atkinson and Sherwood (1992)]; however. 



we did not observe any signature of their existence: they can thus be considered as second-order 
phenomena. 
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Figure 5: Radial velocity profile Vr{r, 6) of a Newtonian oil sheared at 1 rpm in a six-bladed vane- 
in-cup geometry, for various angular positions 9 between two adjacent blades (from ^ = ° to 
^ = 60°). 



As explained in Sec. II, from the Ve(r, ^) measurement and from the material incompress- 
ibility, we are able to reconstruct the radial velocity profile K-(r, 9) (see Fig. 5). This also allows 
us to compute the strain rate components dre{r,9) and drr{r,6) = —dQg{r,6), which are plotted 
in Fig. 6. Of course, due to the limited number of profiles between two adjacent blades, this 
method provides only a rough estimate of these quantities. In addition to their interest for fu- 
ture comparison with models and simulations, these data allow us to evaluate the contribution 
of the extensional flow to dissipation; here, in a Newtonian medium, the local power density is 
given by: pd{r, 9) = 2r]{d% + 2d^^.). 
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Figure 6: Strain rate profiles dr0{r,6) (left) and drrir,9) (right) vs. radius r for various angular 
positions 9 between two adjacent blades (from ^ = ° to ^ = 60 ° ). 



In Fig. 5, we first observe that Vr{r,9) ~ for 9 = ° and 9 = 30°; there is thus no 
extensional flow in these regions of space, as seen in Fig. 6. This is actually expected from the 
fore-aft symmetry of the flow around these angular positions. Vr and its spatial variations (i.e. 
drr) are maximal at ~ 15 ° . Meanwhile, we observe that dre is maximal near the blades: at 
r ~ i?j it is more than 4 times larger at 9 = 0° than at = 30 ° . We then find that dro (and thus 
the shear stress r^g) decreases more rapidly from the blades (at ^ = ° ) than the 1/r^ scaling of 



the Couette geometry, whereas it does not vary much with r midway between adjacent blades 
(it even seems to shghtly increase with r as already observed in simulations by Savarmand et al. 
(2007)). It is also worth nothing that at r ~ i^e; in contrast with what is observed at r c^ Ri, 
the shear stress value is of order two times lower at ^ = ° than at ^ = 30 ° . 

From the whole set of d^e a-nd drr measurements (Fig. 6) , we finally find that in regions where 
drr is maximal, the contribution of the extensional flow to dissipation is of order 25%. Over the 
whole gap, we then evaluate its average contribution to dissipation to be rather important, of 
order 5 to 10%. This significant value may be a reason why the torque that has to be exerted 
to enforce flow is higher than that predicted by Atkinson and Sherwood (1992) in an infinite 
medium. The confinement effect induced by a close boundary at a radius Re likely increases the 
contribution of the extensional flow to dissipation as compared to the case of an infinite medium 
(although other effects may exist, as appears from the comparison of the data of Tab. 1). 

B Yield stress fluid 

In this section, we study the flows induced by the vane tool with a yield stress fluid (a concen- 
trated emulsion). We focus on the behavior near the yielding transition, i.e. on low rotational 
velocities 0,. 
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Figure 7: Dimensionless ^-averaged velocity profile Vo{r,Q)/i^Ri of a yield stress fluid (concen- 
trated emulsion) in a six-bladed vane-in-cup geometry for fi ranging from 0.1 to 9 rpm; the 
vertical dashed line shows the radius of the vane; the dotted line is the theoretical profile for a 
rigid body rotation (for r < Ri). 



In Fig. 7, we plot the 0-averaged azimuthal velocity profiles VQ{r) measured at several Q. 
values ranging from 0.1 to 9 rpm, corresponding to macroscopic shear rates varying between 
0.02 and 2 s~^. We first observe that flow is localized: the material is sheared only up to a 
radius Re < Ro- Re is found to increase as Q increases. This is a classical feature of flows of 
yield stress fluids in heterogeneous stress flelds. It has been observed in Couette geometries 
[Coussot (2005); Ovarlez et al. (2008)], where it is attributed to the 1/r^ decrease of the shear 
stress Trd, which passes below Ty at some Rc{^) < Ro at low $7. In this case, when O tends to 0, 
Re tends to Ri and the torque T at the inner cylinder tends to Ty * 2ttR^H. In the vane-in-cup 
geometry, the same argument holds qualitatively thanks to Eq. 8. It implies that the flow has 
to stop inside the gap at low Vt. However, in contrast with the case of the Couette geometry, 
as the whole stress fleld a priori depends on 0., this ^-averaged equation does not provide the 



position of the limit between the sheared and the unsheared material (which will determined at 
the end of this section) . 

We then observe that, although this effect is less pronounced than with a Newtonian material, 
the shear flow still enters the region between the blades, even at the lowest studied Q. Close 
examination of the profiles shows that the material trapped between the blades rotates as a rigid 
body only up to Ri ~ 3.65 cm at fi = 9 rpm, Ri ~ 3.75 cm at i7 = 1 rpm, and Ri ~ 3.85 cm at 
ft = 0.1 rpm. We recall that Ri ~ 3.05 cm with a Newtonian fluid in the same geometry. 
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Figure 8: Azimuthal velocity profile Vg{r,9) of a yield stress fluid (concentrated emulsion) 
sheared at 0.1 rpm in a six-bladed vane-in-cup geometry, for various angular positions between 
two adjacent blades (from ^ = ° to ^ = 60 ° ). The vertical dashed line shows the radius of the 
vane. 



As in Sec. IIIA3, to better characterize the flow fleld, we have performed experiments in 
which we have measured 10 azimuthal proflles between two adjacent blades at 0.1 rpm. In 
Fig. 8, as for a Newtonian fluid, we observe that there is a strong ^-dependence of the velocity 
proflles. The velocity profile that starts near from a blade tip (at 6 = 0°) has a much steeper 
slope than the profile measured midway between adjacent blades (at ^ = 30 ° ); again, this shows 
that the blade tip neighborhoods are regions of high shear. Meanwhile the flow stops at a radius 
Re which is larger at ^ = 30 ° (4.5 cm) than at ^ = ° (4.3 cm). Note also that there may 
be slight fore-aft asymmetry, as sometimes observed with yield stress fluids flows [Dollet and 
Graner (2007); Putz et al. (2008)], but we did not study this point further. From these velocity 
proflles, we have reconstructed a 2D map of the flow field (Fig. 9), indicating both the boundary 
between the region of rigid body rotation (between the blades) and the sheared region, and the 
boundary between the sheared region and the outer region of fiuid at rest (i.e. the position 
where the yield criterion is satisfied). 

Flow is found to occur in a layer of complex shape which is far from being cylindrical even at 
this very low velocity. These observations are in contradiction with the usually accepted picture 
for yield stress fluid flows at low rates [Barnes and Nguyen (2001)], namely that the material 
contained in the virtual cylinder delimited by the blades rotates as a rigid body, and that it 
flows uniformly in a thin cylindrical layer near the blades. Our results contrast in particular 
with previous numerical works which showed that the yield surface is cylindrical at low rates 
for Bingham fluids, Casson fluids, and Herschel-Bulkley materials with n = 0.5 [Keentok et al. 
(1985); Yan and James (1997); Savarmand et al. (2007)]. With apparently similar conditions 
to those in some of the Yan and James (1997) simulations, we find an important departure 
from cylindrical symmetry. This means that further investigation on the exact conditions under 
which this symmetry can be recovered is still needed. Possible difference between our work and 



that of Yan and James (1997) is that the blade thickness is zero in this last study. 

It is particularly striking and counterintuitive that Re is largest at the angular position 
{9 = 30°) where shear at Ri is smallest (similar observation was made by Potanin (2010)). 
As in Sec. IIIA3, this points out the importance of the extensional flow in this geometry, with 
^-dependent normal stress differences which have to be taken into account in the yield criterion, 
and which thus impact the yield surface. It thus seems that the link between the yield stress 
Ty and the torque T measured at yield with a vane-in-cup geometry is still an open question, 
although the classical formula probably provides a sufficiently accurate determination of Ty in 
practice. 
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Figure 9: Two-dimensional plot of the limit between rigid motion and shear (circles) and between 
shear and rest (triangles) for a yield stress fluid (concentrated emulsion) sheared in the six-bladed 
vane-in-cup geometry at 0.1 rpm (left) and 1 rpm (right). The grey rectangles correspond to 
the blades. 



The same 2D map as above is plotted for 0, = 1 rpm in Fig. 9; the same phenomena are 
observed, with enhanced departure from cylindrical symmetry, consistent with the observation 
that Ri decreases when fi increases. This result was also unexpected, as simulations find uniform 
flows for shear-thinning material of index n < 0.5 [Barnes and Carnali (1990); Savarmand et 
al. (2007)]; we would have expected the same phenomenology in a Herschel-Bulkley material of 
index n = 0.5 (and thus Ri to tend to Ri when increasing Vi). This observation also shows that 
a Couette analogy can hardly be defined for studying the flow properties of such materials in 
a vane-in-cup geometry because the equivalent Couette geometry radius Ri^cq would probably 
depend also on 17 (as recently shown by Zhu et al. (2010)). 

Let us finally note that this departure from cylindrical symmetry has important impact on 
the migration of particles in a yield stress fluid (see below) . 



C Concentrated suspension 

In this section, we investigate the behavior of a concentrated suspension of noncolloidal particles 
in a yield stress fluid (at a 40% volume fraction). 

A detailed study of their velocity profiles would a priori present here limited interest: such 
materials present the same nonlinear macroscopic behavior as the interstitial yield stress fluid, 
and their rheological properties (yield stress, consistency) depend moderately on the particle 
volume fraction [Mahaut et al. (2008a); Chateau et al. (2008)]. 

On the other hand, noncolloidal particles in suspensions are known to be prone to shear- 
induced migration, which leads to volume fraction heterogeneities. This phenomenon is well 



documented in the case of suspensions in Newtonian fluids [Leighton and Acrivos (1987b); 
Abbott et al. (1991); Phillips et al. (1992); Corbett et al. (1995); Shapley et al. (2004); Ovarlez 
et al. (2006)] but is still badly known in yield stress fluids (some studies exist however in 
viscoelastic fluids [Tehrani (1996); Huang and Joseph (2000); Lormand and Phillips (2004)]). 
In the model of Leighton and Acrivos (1987b) and Phillips et al. (1992), migration is related 
to shear-induced diffusion of the particles [Leighton and Acrivos (1987a); Acrivos (1995)]. In a 
wide gap Couette geometry, the shear stress heterogeneity is important (Eq. 4); the shear rate 
gradients then generate a particle flux towards the low shear zones, which is counterbalanced 
by a particle flux due to viscosity gradients. A steady state, which results from competition 
between both fluxes, may then be reached, and is characterized by an excess of particles in the 
low shear zones of the flow geometry (near the outer cylinder in a wide-gap Couette geometry 
[Phillips et al. (1992); Corbett et al. (1995); Ovarlez et al. (2006)]). Note that there are other 
models [Nott and Brady (1994); Mills and Snabre (1995); Morris and Boulay (1999); Lhuilher 
(2009)] in which particle fluxes counterbalance the gradients in the particle normal stresses, and 
which can be used directly for non-Newtonian media. 

As the development of migration depends on the spatial variations of shear, one may wonder 
how the azimuthal heterogeneities of shear introduced by the vane tool affect migration; a related 
question is that of the relevance of the Couette analogy for this phenomenon. In the following, 
we thus focus on the particle volume fraction distribution evolution when the material is sheared. 

Behavior at high shear rate 

We first study the behavior at high shear rate, in the absence of shear localization. We shear the 
suspension in both the standard Couette geometry and the vane-in-cup geometry at a rotational 
velocity J7 = 100 rpm. In this first set of experiments, we only study the steady-state of 
migration. At Jl = 100 rpm, this steady-state is reached in less than 30 min (which corresponds 
to a macroscopic strain of order 50000, consistently with strainscale evaluations from data of the 
literature [Ovarlez et al. (2006)]). In Fig. 10 we plot the steady state volume fraction profiles 
observed after shearing the suspension at i7 = 100 rpm during Ih. 




(a) 



3 4 

Radius (cm) 




(b) 



3 4 5 6 

Radius (cm) 



Figure 10: a) Steady-state volume fraction vs. radius at Q = 100 rpm in both the Couette 
geometry (empty circles) and the vane-in-cup geometry (squares) . In the vane-in-cup geometry, 
the volume fraction profile is determined in a 1 cm thick slice situated exactly between two 
adjacent blades (see Fig. lib). The inset is a zoom; the line is a fit of the data measured in 
the Couette geometry to the Phillips et al. (1992) model with Kc/K^ = 0.42. b) 0-averaged 
azimuthal velocity profile Ve(r); the dotted line is the theoretical rigid motion induced by the 
rotation of the vane tool; the vertical dashed line shows the radius of the vane. 



As expected, we first observe that the material is strongly heterogeneous in the Couette 
geometry: the volume fraction varies between 37% near the inner cylinder and 43% near the 
outer cylinder (note that the NMR technique we use do not allow quantitative measurements 
near the walls). This heterogeneity is quantitatively similar to that observed in Couette flows of 
Newtonian suspensions at a same 40% particle volume fraction [Corbett et al. (1995)]; the profiles 
are actually well fitted to the Phillips et al. (1992) model (see Eq. 16 of Ovarlez et al. (2006)) 
with a dimensionless diffusion constant Kc/K^ = 0.42 which is close to that found by Corbett 
et al. (1995) {Kc/K^ = 0.36), although this model is not expected to hold in non- Newtonian 
suspensions. 

In the vane-in-cup geometry, the volume fraction profile shows very different features; note 
that the profile is measured in a 1 cm thick (in the azimuthal direction) slice situated exactly 
between two adjacent blades (see Fig. 11). In Fig. 10, we first observe that there is a strong 
particle depletion in a wide zone between the blades. A homogeneous volume fraction of 40% 
is observed for radii inferior to 3.1 cm. At a radius Ri = 3.1 cm, there is a strong drop in 
the volume fraction down to 5% within 1 mm (corresponding to 4 particle diameters). Close 
inspection of the velocity profile Fig. 10 shows that this radius Ri corresponds to the transition 
between rigid motion and shear between the blades. The volume fraction then increases basically 
linearly with the radius up to a 40.5% volume fraction at a radius r = 3.85 cm which is close to 
the vane radius. The volume fraction then increases only slightly (between 40.5% and 42.5%) 
in the gap of the geometry: the heterogeneity is here much less important than in a standard 
Couette geometry. 

To get further insight into the new strong depletion phenomenon we have evidenced, we have 
performed 2D magnetic resonance images of the material. Such images provide a qualitative 
view of the spatial variations of the particle volume fraction as only the liquid phase is imaged. 
Images are coded in grey scales; a brighter zone contains less particles. In Fig. lib, we first 
see an image of the homogeneous material. Before any shear, as expected, the light intensity 
is homogeneous in the sample (intensity variations correspond to noise). After a Ih shear at 
il = 100 rpm, we observe very bright and thin curves on the image: they correspond to zones 
where the volume fraction suddenly drops down to a value close to zero. These curves are not 
circles. More precisely, between two adjacent blades, a depleted zone goes from the edge of one 
blade (at 6* = ° , r = 4.02 cm) to the edge of another blade (at = 60 ° , r = 4.02 cm), and 
describes a concave r{9) curve whose minimum is r = 3.1 cm at = 30 ° . Note that as the 
volume fraction profile is averaged over a slice which is 1 cm thick in the azimuthal direction 
(see Fig. lib), the fact that we measure a minimum of 5% at r = 3.1 cm in the slice probably 
means that the volume fraction minimum is actually equal to zero in the depletion zone. 

As pointed out above, this curve also likely marks the transition between the unsheared 
material (which rotates as a rigid body) and the sheared material. Note in particular the 
similarity with Fig. 4b, the data of which are reported in Fig. 11a for illustration. A first 
interpretation of the phenomenon would then simply be that migration is caused by shear and 
naturally stops at this transition zone. Indeed, as shear is maximum near the blades, particles 
tends to migrate out of this zone; moreover, there is no source of particle fiux from the unsheared 
zone between the blades to balance the migration towards the outer cylinder. However, this does 
not explain why the volume fraction drops down to zero: heterogeneities observed at steady- 
state in the literature are usually moderate and do not lead to zones free of particles. A better 
understanding of the phenomenon can be gained by zooming on the previous image (Fig. lie). 
We now see that while the depletion phenomenon seems symmetric around both sides of the 
blades at a macroscopic scale, it is clearly asymmetric at a local scale near the blades and depends 
on the direction of rotation: depletion is more pronounced at the back of the blade (note that the 
vane tool rotates counterclockwise). This would mean that the noncolloidal particle trajectories 
are asymmetric around the blade: a particle that is found at a radius r ~ 4.02 cm just before 
the blade is necessarily found at a radius slightly higher than 4.02 cm after the blade as there 




(a) 




Figure 11: 2D magnetic resonance image of a suspension of particles in a yield stress fluid in 
a vane-in-cup geometry: (a) after a Ih shear at fi = 100 rpm (corresponding to a macroscopic 
strain of order 75000), and (b) before any shear. The crosses in Fig. 11a correspond to the limit 
between rigid motion and shear for the Newtonian oil of Fig. 4b. The white rectangle in Fig. lib 
shows the slice in which the volume fraction profiles of Figs. 10a and 12 are measured, (c) is a 
zoom of image (a) near the edges of a blade. The images are taken in the horizontal plane of 
the geometry, at middle height of the vane tool, and correspond to vertical averages over 2 cm. 
The vane tool rotates counterclockwise. 



are no particles at r = 4.02 cm. This feature is reminiscent of the fore-aft asymmetry that is 
observed in the bulk of noncolloidal suspensions [Parsi and Gadala-Maria (1987)] and that leads 
to their non-Newtonian properties [Brady and Morris (1997)]. It thus seems that, in addition to 
the shear-induced migration mechanism intrinsic to suspensions, the vane tool induces a specific 
migration mechanism which has its origin in the direct interactions between the particles and 
the blades; this effects leads to the full depletion that is observed at the transition between the 
sheared and the unsheared material. Such direct effect of a flow geometry on migration has also 
been observed in microchannel flows of colloidal suspensions [Wyss et al. (2006)], and also led 
to full particle depletion. See also Jossic and Magnin (2004). The kinetics of the phenomenon 
will be briefly discussed below. 

The rest of the volume fraction profile results from a complex interplay between shear- induced 
migration and the fore-aft asymmetry around the blades; this leads to the rapid increase of the 
volume fraction between 3.1 cm and 4.02 cm. After 4.02 cm the fiow lines do not meet the blade 
edges, and the phenomenon evidenced above should have basically no effect on the heterogeneity 
that develops in the gap of the geometry. On the other hand, the mean volume fraction should 
be slightly higher due to mass conservation; it is indeed observed to be equal to around 42%. 
Nevertheless, as the mean radial shear rate heterogeneity is basically similar to that observed in 
a standard Couette geometry (see previous sections), we would a priori expect the heterogeneity 
to be somehow similar. However, we observe that the volume fraction profile is only slightly 
heterogeneous: there is less than 5% variation of the volume fraction in the gap, to be compared 
to the 15% variation observed in the Couette geometry. Clearly, this means that the Couette 
analogy is irrelevant as regards this phenomenon, and that the details of shear matter. Here, 
the extensional flow that adds to shear may be at the origin of this diminution of migration. A 
more detailed analysis is out of the scope of this paper. 

Behavior at low shear rate 

Let us now study the behavior at low shear rate. Low shear rates are typically imposed with the 
aim of measuring the yield stress of such materials. Starting from a homogeneous suspension 
at rest, we apply a rotational velocity = 1 rpm (without any preshear), and we measure the 
evolution of the particle volume fraction in time. The corresponding volume fraction proflles are 
depicted in Fig. 12. 

In Fig. 12, we observe that, although shear is much less important than in the previous 
experiments, particle depletion also appears between the blades. Comparison of the velocity 
profile and the volume fraction profile shows that depletion also appears between the blades 
at the transition zone between the sheared and the unsheared materials. This phenomenon 
appears with a very fast kinetics: the lower volume fraction value in the measurement zone is 
36% after only a 5 minute shear (corresponding to a macroscopic strain of order 50). Afterwards, 
it continues evolving slowly: the minimum observed volume fraction is of order 33% after a Ih 
shear and of order 32% after a 14h shear (corresponding to a 10000 macroscopic strain). Note 
that the radial position of the minimum value of the volume fraction slightly decreases in time; it 
likely corresponds to progressive erosion of the material between the blades (we did not measure 
the velocity profiles to check this hypothesis). 

We also note that migration is negligible in the rest of the sheared material as expected 
from the theory of migration briefly described above (a larger strain would be needed to observe 
significant migration). Nevertheless, we note some particle accumulation (with a volume fraction 
value of 43%) at i?c=4.7 cm after a very long time. This corresponds to the yield surface as flow 
is localized at low velocity (see velocity proflle Fig. 12). Migration profiles usually result from an 
equilibrium between various sources of fluxes. On the other hand, the unsheared material does 
not produce any particle flux while it receives particles from the sheared region. This particle 
accumulation is thus the signature that the migration phenomenon is indeed active, although 
not observable on the profile measured in the sheared zone. It is probable that this accumulation 
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Figure 12: Volume fraction vs. radius at il = 1 rpm measured in the vane-in-cup geometry after 
different times of shear: 5min, Ih, 14h. The material is homogeneous at the beginning of shear. 
The inset presents the ^-averaged azimuthal velocity profile Voir) measured in the first stages 
of shear; the dotted line is the theoretical rigid motion induced by the rotation of the vane tool; 
the vertical dashed line shows the radius of the vane. 



process would stop only (after a very long time) when there are no more particles in the sheared 
region. 

As above, 2D magnetic resonance images of the material provide an insight in the phe- 
nomenon. In Fig. 13, we observe again that particle depletion is enhanced at the rear of the 
blades; this confirms that this phenomenon is likely due to direct interactions between the blades 
and the particles, leading to the asymmetry of the particles trajectory around the blades. This 
a priori occurs with any particle whose trajectory is close to the blades, explaining why particle 
depletion appears so rapidly. There is probably no way to avoid it. Note that the images are 
here much brighter very close to the blades than midway between two adjacent blades; this 
would mean that the particle volume fraction is probably close to near the blades, although 
we observe volume fraction of order 32% between two blades. 

Finally, let us note that the bright line provides a good idea of the boundary between the 
sheared material and the material that moves as a rigid body. We see as in Sec. IIIB that this 
is far from being cylindrical even at this low velocity. 



Consequences: slip ^vith a vane tool 

We finally present some consequences of this phenomenon. From the above observations, our 
conclusion is that depletion sets up quickly and is probably unavoidable. Then two situations 
have to be distinguished. If linear viscoelastic properties of a suspension of large particles are 
measured at rest on the homogeneous material in its solid regime, without any preshear, then 
these measurements pose no other problem than that of the relevant Couette analogy to be used 
(see Sec. IIIA). If a yield stress measurement is performed at low imposed rotational velocity, 
starting from the homogeneous material at rest, then this measurement is likely valid as long as 
only the peak value or the plateau value at low strain (of order 1) is recorded. On the other hand, 
any subsequent analysis of the material behavior will a priori be misleading: irreversible changes 
have occurred and the material cannot be studied anymore. More generally, any measurement 
performed after a preshear will be incorrect and any flow curve measurement will lead to wrong 
evaluation of the material properties. In these last cases, the consequence of the new particle 
depletion phenomenon we have evidenced is a kind of wall slip near the blades, whereas there 




Figure 13: (a) 2D magnetic resonance image of a suspension of particles in a yield stress fluid 
in the horizontal plane of a vane-in-cup geometry after a 14h shear at = 1 rpm. (b) is a zoom 
of image (a) between two adjacent blades. The image is taken in the horizontal plane of the 
geometry, at middle height of the vane tool, and corresponds to a vertical average over 2 cm. 
The vane tool rotates counterclockwise. 



are no walls. Here the "slip layer" is made of the (pure) interstitial yield stress fluid in a zone 
close to the blades, as would be observed near a smooth inner cylinder. This contrasts with the 
common belief that the vane tool prevents slippage. 

In order to illustrate this feature, we present some results of Mahaut et al. (2008a): Mahaut 
et al. performed classical upward/downward shear rate sweeps with a six-bladed vane-in-cup 
geometry [Ri = 1.25 cm, Rq = 1.8 cm, H =4.5 cm, blade thickness=0.8 mm) in a pure con- 
centrated emulsion, and in the same emulsion filled with 20% of 140 jiva. PS beads. In these 
experiments, constant macroscopic shear rates steps increasing from 0.01 to 10 s~^ and then 
decreasing from 10 to 0.01 s~^ were applied during 30s, and the stationary shear stress was 
measured for each shear rate value. The results are shown in Fig. 14. 

While the same curve is observed for the upward/downward shear rate sweeps in the case of 
the pure emulsion (as expected for a simple non-thixotropic yield stress fluid), the shear stress 
during the upward shear rate sweep differs from the shear stress during the downward shear rate 
sweep in the case of the suspension. Moreover, any measurement performed on the suspension 
after this experiment gives a static yield stress equal to the dynamic yield stress observed during 
the downward sweep. This means that there has been some irreversible change. This irreversible 
change is actually the particle depletion near the blades we have observed in this paper. As the 
flow of the suspension is localized near the inner tool at low shear rate, it means that after 
the first upward sweep that has induced the particle depletion, during the downward shear rate 
sweep only the pure emulsion created by migration near the blades remains in the sheared layer 
at sufficiently low rotational velocity. This explains why the same apparent value of the yield 
stress is found in the suspension during the downward sweep as in the pure emulsion with this 
experiment. On the other hand, the yield stress at the beginning of the very first upward sweep 
is that of the suspension as migration has not occurred yet. 

The conclusion is that the vane tool is probably not suitable to the study of flows of suspen- 
sions of large particles. 
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Figure 14: Shear stress vs. shear rate in a vane-in-cup geometry for upward/downward shear 
rate sweeps in a pure concentrated emulsion (open squares) and for the same emulsion filled 
with 20% of 140 fim PS beads (filled/open circles). Figure from Mahaut et al. (2008a). 



IV Conclusion 

As a conclusion, let us summarize our main findings: 

• In the case of Newtonian fluid flows, our measurements support the Couette equivalence 
approach: the 0-averaged strain rate component d^e decreases as the inverse squared radius 
in the gap. Interestingly, the velocity profiles allow determining the Couette equivalent 
radius without end-effect correction and independently of the viscosity of the material. 
The torque exerted by the vane in our display is found to be higher (by 8%) than the 
theoretical prediction of Atkinson and Sherwood (1992) for a vane embedded in an infinite 
medium, and is thus much closer to the torque exerted by a Couette geometry of same 
radius as the vane than expected from the literature. A key observation may be that there 
is a significant flow between the blades which adds an important extensional component to 
shear, thus increasing dissipation. From a short review of the literature, it clearly appears 
that numerical investigations are still needed in the case of finite geometries. Variational 
approaches are also promising, although they do not yet provide tight bounds. 

• In the case of yield stress fluid flows, we find that the thin layer of material which flows 
around the vane tool at low velocity is not cylindrical, in contrast with what is usually 
supposed in the literature from simulation results. Consequently, a non negligible exten- 
sional component of shear has probably to be taken into account in the analysis. At this 
stage, there are too few experimental and simulation data to understand the origin of this 
discrepancy. It thus seems that progress still has to be made, in particular through simula- 
tions, which allow a wide range of parameters to be studied. This may help understanding 
how the torque is linked to the yield stress of a material at low velocity, depending in 
particular on the geometry. 

• An important and surprising result is the observation of particle depletion near the blades 
when the yield stress fluid contains noncolloidal particles. This phenomenon is thus likely 
to occur when studying polydisperse pastes like coal slurries, mortars and fresh concrete. 
It has to be noted that the phenomenon is very rapid, irreversible, and thus probably 
unavoidable when studying flows of suspensions of large particles. It results in the creation 
of a pure interstitial yield stress fluid layer and thus in a kind of wall slip near the blades. It 
contrasts with the classical assumption that is made in the field of concentrated suspension 
rheology where the vane tool is mainly used to avoid this phenomenon. 



Consequently, we would say that, in the case of pasty materials with large particles, if accurate 
measurements are needed, the vane tool may finally be suitable only for the study of the solid 
(elastic) properties of materials and for the static yield stress measurements; as the yield stress 
measurement may induce irreversible particle depletion near the blades, any new measurement 
then requires a new sample preparation. Furthermore, the vane can be used as a very accurate 
tool without any hypothesis nor any calibration to measure the relative increase of the elastic 
modulus of materials as a function of their composition [Alderman et al. (1991); Mahaut et al. 
(2008a)]. In order to study accurately the flows of pasty materials with large particles, our results 
suggest that a coaxial cylinders geometry with properly roughened surfaces is preferable when 
possible. If the use of a vane tool cannot be avoided, one should keep in mind our observations 
in order to carefully interpret any result. 
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